Edge States and Quantum Hall Effect in Graphene under a Modulated Magnetic Field 
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Graphene properties can be manipulated by a periodic potential. Based on the tight-binding 
model, we study graphene under a one-dimensional (ID) modulated magnetic field which contains 
both a uniform and a staggered component. The chiral current-carrying edge states generated at 
the interfaces where the staggered component changes direction, lead to an unusual integer quantum 
Hall effect (QHE) in graphene, which can be observed experimentally by a standard four-terminal 
Hall measurement. When Zeeman spin splitting is considered, a novel state is predicted where the 
electron edge currents with opposite polarization propagate in the opposite directions at one sample 
boundary, whereas propagate in the same directions at the other sample boundary. 

PACS numbers: 71.70.Di, 73.43.Cd, 73.61.Wp 



Recently, graphene materials have received extensive 
theoretical and experimental studies The most im- 
portant physical properties of graphene are governed by 
the underlying chiral Dirac fermionsQ, These Dirac 
fcrmions under a uniform magnetic field (UMF) give 
rise to the well-known anomalous QHE0, Hj], which has 
been experimentally verified and is now believed to be 
a unique feature to characterize graphene. Spin QHE 
was also predicted in graphene, where electron edge cur- 
rent with the opposite spin polarization couterpropa- 
gates due to the spin-orbit interaction Q or the Zee- 
man spin splitting of the zcroth Landau level (LL)0j. 
On the other hand, the experimental manipulation of 
the electronic structure of graphene has potential ap- 
plication in graphene electronics or spintronicsj^, 
One method to manipulate the physical properties of 
graphene is by applying periodic electronic [3, [ll[ or 
magnetic potentials Il3j . which can be realized now 
by making use of substrate [3 M, M, E3 

or controlled 

adatom deposition [ijj]. 

Here, we report the investigation of the effect on 
graphene QHE of a ID staggered magnetic field (SMF), 
which is schematically shown in the top panels of Fig. [1] 
The ID SMF can be achieved in experiments by applying 
an array of ferromagnetic stripes with alternative mag- 
netization on the top of a graphene layer. It is found 
that the edge states created by ID SMF lead to a non- 
trivial robust integer QHE in graphene. In a standard 
four-terminal Hall measurement, when varying the mag- 
nitude of the UMF, graphene can undergo a transition 
from a state with unusual quantized Hall conductance to 
one without Hall effect. Furthermore, the Zeeman spin 
splitting of the zeroth LL of graphene gives rise to a novel 
state where spin-up and spin-down edge currents have the 
opposite chirality at one sample boundary but have the 
same chirality at the other sample boundary. 

We start from the tight-binding model on a honeycomb 
lattice in a perpendicular nonuniform magnetic field dc- 
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FIG. 1: (color online). Top panels: Schematic illustration of 
a rectangular graphene sample under a SMF, where two uni- 
form fields which are the same in magnitudes but opposite 
in directions alternate every Lx distance. Case I in (a) con- 
tains 2N X alternating regions whereas Case II in (b) contains 
2N X + 1 alternating regions, with N x an integer. Bottom pan- 
els: Electron energy bands of graphene under a periodically 
SMF with L x = 426 nm, and B s = 24 T for (c), B s = 40 
T for (d). Energy is measured in unit of KvfIIb with Ib the 
magnetic length. For the parameters chosen, Ib is determined 
to be Ib = 3.7 nm for (c), Ib = 2.9 nm for (d), respectively. 
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where t is the hopping integral, the operator cj (a) cre- 
ates (annihilates) an electron at site i, and < ij > de- 
notes nearest-neighbor pairs of sites. A is the gauge po- 
tential for the nonuniform magnetic field. The Zeeman 
spin splitting is neglected now for simplicity and will be 
discussed later. We distinguish here Case I in Fig. [TJa) 
from Case II in Fig. |Tfb) because graphene in Case I has 
no QHE unless the magnetic fields of the two alternat- 
ing areas have the same directions, i.e. Bjj > Bs where 
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Bu and B$ are the magnitudes of the UMF and SMF 
respectively. The reason for this will be explained later 
when we discuss Fig. [2] 

The LLs of graphene [lj| can be expressed as E n = 
±HvF\/\n\/lB with Ib = \Jhcj2eB the magnetic length, 
vp = 3at/2h the Fermi velocity, and n = 0, 1, 2, ... the LL 
index. Here a is graphene lattice constant. The physical 
picture at large L x order of L x >> Is is found to be 
quite different from that at smaller L x case with L x ~ Ib, 
where only a few chains are contained in each alternating 
area. In the latter case, though Dirac cone structure is 
preserved, more and more Dirac points are created as 
increasing L^(20|. and finally the LLs of graphene appear. 
In Figs.QJc) andQJd) we show in the absence of the UMF 
the electron energy bands of graphene under a periodic 
SMF with k x = 0. For the SMF with a period 2L X , we 
have chosen correspondingly a periodic gauge A = A# = 
(0,Bs\x\,0) for \x\ < L x . When energy is set by Kvf/Ib, 
the \fn spacing of the LLs of graphene can be seen clearly. 
What's remarkable is that for large L x , the dense Dirac 
points are emerged into the zeroth LL of graphene (20j. 
Compared with the LLs of graphene in a UMF, which 
is dispersionless, the energy bands have non-flat regions 
where energy disperse with k y , indicating the presence 
of edge states. These new edge states which were called 
magnetic edge states 21, [22| are actually generated by 
the SMF and right located at the interfaces where the 
SMF changes direction, providing the gapless excitations 
for even an infinite graphene. 

To further clarify the physical properties of these edge 
states and its consequences, we study the QHE of such a 
system, i.e., the graphene samples in the presence of both 
a UMF and a SMF. In the following cases, open bound- 
ary condition is applied in the x direction and periodic 
boundary condition in the y direction, and the Landau 
gauge Ay = (0, Bjjx, 0) is adopted for the UMF and 
A = As + Ay. The graphene energy spectrum for Case 
I and Case II are calculated and shown in Figs. [2ja) 
and[2]Jb), respectively. With application of a UMF, the 
graphene samples can be divided into two groups of re- 
gions, where one is the high-field group of regions with 
the magnetic field Bjj + Bs, the other is the low-field 
group of regions with the magnetic field Bu — B$- Cor- 
respondingly, the bulk excitations of the graphene sys- 
tem have two groups of LLs, which are reflected in elec- 
tron energy bands by the two groups of the flat regions 
in Figs. [2ja) and [UJb). Interestingly, when the ratio 
(Bu — Bs)/(Bu + Bs) = p/q, where p and q are two 
coprime integers, a series of LLs will be doubly degen- 
erate, which may cause some interesting phenomena. In 
particular, we actually have (Bu — Bs) / (Bu + Bs) = 1/5 
and —1/3 for the parameters in Figs.^a) and[2jb), re- 
spectively, resulting in the doubly degeneracy of all the 
LLs in the high- field regions. 

Now we focus our attention on the edge states. In 
Figs. (2)Jc) and [U^d), electron probability densities for 
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FIG. 2: (color online) . Top panels: Electron energy spectrum 
of graphene in the presence of both a UMF and a SMF for 
N x — 1 and L x = 213 nm. (a) belongs to Case I with pa- 
rameters Bu = 24 T, B s = 16 T. (b) belongs to Case II with 
parameters Bu = 16 T, Bs — 32 T. Here energy is scaled 
in the same way to Fig. [1] except the magnetic length / b is 
replaced by that for the high-field regions. The red dashed 
lines represent the positions of the chemical potentials. Mid- 
dle panels: The corresponding electron probability densities 
for the representative states indicated in the top panels by 
the open circles. Bottom panels: Schematic diagrams of the 
corresponding electron edge currents (green and blue arrows) 
and electrons' classical skipping orbits. Vh is the Hall voltage. 
Gray (light gray, yellow) solid rectangle represents high- (low-) 
field regions. Black and white solid circles represent particle- 
hole excitations in the Laughlin's thought experiment. 



representative edge states are shown. Except the con- 
ventional edge states located at the sample boundaries 
[see (1),(2),(4),(6) in Fig. He), and (1),(2),(7),(8) in 
Fig. Hfd)], new edge states are generated and right lo- 
calized at the interfaces where the SMF changes direc- 
tion [see (3), (5) in Fig. Mc), and (3)-(6) in Fig. E^d)]. 
Clearlvjthese edge states are also chiral current-carrying 
states [H, [U [25[, whose flowing directions can be easily 
determined by the slopes of the bands where they are 
located. All the edge currents are schematically shown 
in the bottom panels of Fig. [2] Also indicated in these 
panels are the classical orbits of electrons}^, which arc 
composed of arcs with the length scale approximately 
equal to the magnetic length Is- These classical orbits 
give a simple physical interpretation of the reason why 
the edge currents flow in the shown directions. 

Let us now consider the corresponding four-terminal 
Hall measurements. Assuming that contacts are reflec- 
tionlcss, all edge currents coming from the same contact 
share the same voltage with that contact, so that all the 
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FIG. 3: (color online). The four-terminal Hall bar geometries 
for (a) Case I and (b) Case II with N x = 3. The black 
solid rectangles represent contacts and the lines with arrows 
represent the edge currents. 



"up-flowing" edge currents have the same voltage and 
so do all the " down- flowing " edge currents. The volt- 
age difference between the edge currents with opposite 
directions is equivalent to the Hall voltage Vh between 
the sample boundaries. To obtain the Hall conductance 
in these two situations in Figs.[^e)-(f), we follow Laugh- 
lin's argument (27j and imagine rolling up the graphene 
ribbon along the y direction to make a graphene cylinder 
which is then threaded by a magnetic flux. When vary- 
ing adiabatically the magnetic flux through the graphene 
cylinder by a flux quantum </>o, particle-hole excitations 
will be generated both at the sample boundaries and the 
interfaces where the SMF changes direction. For an il- 
lustration, we show schematically in Figs. Ele)-(f) these 
excitations for the chemical potentials indicated by the 
red lines in Figs. [2]^a)-(b). The increase in energy due to 
electrons transfer between the Hall voltage is SE = 3eVn 
for both Figs.[2jc) and[2^f). Generally, if the system has 
such a Fermi energy that the Mth LL of the high-field 
regions and the Nth LL of the Low-field regions are just 
completely filled, detailed analysis of the edge states in 
the energy bands shows that the increase in energy SE 
can be given by SE = P^eVa for Case I in Fig.^e), and 
SE = (P N + 2P M )eV H for Case II in Fig. [^f), respec- 
tively, where P M = 2M + 1 and P N = 2N + 1. Hence, 
from I = cSE/Stfi, the Hall conductance is given by 
<?xy = Pnc 2 /h for Case I in Fig.^c), which has nothing 
to do with the high-field LLs, and a xy = (Pjy + 2PM)e 2 /h 
for Case II in Fig.^f). 

We note that in order to have the mentioned QHE for 
Case 1(H), the magnetic field for the high-field and low- 
field regions must have the same (opposite) direction, i.e., 
Bu > B s {Bjj < B s ). For Case I in Fig.^e), by varying 
the UMF so that the magnetic field for the low-field re- 
gion is reversed, i.e., Bu < Bs, the electron edge currents 
at the right edge and the interface will be reversed too, 
resulting in the same voltage shared by the two sample 
boundaries and thus a zero Hall voltage. For Case II in 
Fig- [SJf ) , however, the reversal of the magnetic field in 
the low-field region for Bjj > Bs only leads to the rever- 
sal of edge currents at the two interfaces with that at the 
two sample boundaries unchanged, giving rise to another 
quantized Hall conductance a xy — P^e 2 /h. Therefore, 



for both cases, Bu — Bs is a critical value for graphene 
QHE under a SMF. Another important point we should 
remark is that for Case II in Fig. [2c), even in the absence 
of a UMF, graphene under a purely SMF shows a quan- 
tized Hall conductance a xy = ZPjqe 2 jh, where P^ = Pm- 
This can be comparable with the Haldane model intro- 
duced in Ref . [28| , where there exists a nonzero quantized 
Hall conductance in the absence of an external UMF. All 
these peculiar behaviors are believed to have great appli- 
cation in graphene manipulation. 

It is straightforward to generalize the scheme used 
above to more general and interesting cases for N x > 1. 
Our detailed calculation confirm the existence of the two 
groups of the LLs for both Cases, which are the bulk ex- 
citations of graphene under the modulated magnetic field 
and are represented by the flat bands in the electron spec- 
trum. Except the LLs and conventional edge states at 
the sample boundaries, there exist many current-carrying 
edge states localized at the interfaces where the SMF 
changes direction. In Fig. [3l we take N x = 3 for ex- 
ample, and show schematically the corresponding four- 
terminal measurements, where N x now has the meaning 
of the number of the low-field regions. Detailed analy- 
sis of these edge states from the spectrum and similar 
argument lead to the Hall conductance as follows: 
when Bjj > Bs, for both Case I and II, 

a\ y = all = ±j[N x P N -(N x -l)P M ], (2) 
when Bu < Bs, for Case II, 

a x \ = ±jW*P*r + (N x + 1)Pm], (3) 

whereas for Case I, there is no Hall voltage. Here the 
"±" symbol represents the particle-hole symmetry. If 
N x = 1, the previous results are recovered. The result 
Eq. (|3|) seems trivial since it can be seen as the con- 
ductance sum of all the contributions from each distinct 
region. The result Eq. ©, however, is highly nontrivial 
since it can not be seen as the naive subtraction of the 
contribution of the high-field regions from that of the 
low-field regions. When the chemical potential is within 
a gap between the LLs, while the Hall conductance is 
quantized, the resistance of system should be vanishingly 
small because only the edge states carry current so that 
the backscattering is strongly suppressed. Due to the 
gauge invariance of Laughlin's argument used here, the 
results we obtained should be robust against weak disor- 
der and interaction. On the other hand, by varying the 
magnitude of the UMF, we come to a similar conclusion 
to the cases with N x = 1 that, the graphene will undergo 
a transition, from one state with the quantized Hall con- 
ductance Eq. ((2]) to one without Hall effect for Case I in 
Fig. EJa), whereas from one state with the Hall conduc- 
tance Eq. © to one with the Hall conductance Eq. @ 
for Case II in Fig.^b). 
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FIG. 4: (color online). Top panels: Electron energy spectrum 
of graphene in the presence of Zeeman splitting for (a) Case 
I and (b) Case II with N x = 1. The Zeeman energy is cho- 
sen to be one tenth of the largest LL spacing of the high-field 
regions. The red dashed lines represent the positions of the 
chemical potentials. The blue filled circles indicate the edge 
states at the sample boundaries whereas the green filled cir- 
cles indicate the newly generated edge states localized at the 
interfaces where the SMF changes direction. Bottom panels: 
The corresponding four-terminal Hall bar geometries, where 
the red and black lines with arrows represent the spin-up and 
spin-down edge currents, respectively. 



Now we turn to explore the spin effect in graphene 
QHE in the presence of Zeeman spin splitting. The 
energy spectrum is shown in Figs. H[a)-(b). Two spin 
gaps Ah and A^, which are corresponding to the Zee- 
man splitting in the high-field and low-field regions re- 
spectively, are opened. When the chemical potential lies 
within the interval — A^/2 < /j, < Al/2, there is no edge 
current at the interfaces where the SMF changes direc- 
tion, as well as two edge states at both sample bound- 
aries, which have opposite spin polarizations and propa- 
gate in opposite directions [7[. 

When the chemical potential lies within the interval 
Al/2 < \fj,\ < Ah/2, fully spin-polarized edge currents 
appear at the interfaces where the SMF changes direc- 
tion (See the four-terminal Hall measurements shown 
in Figs. IUc)-(d)), which is spin-down for Bjj > Bs, 
and spin- up for Bjj < Bs- Remarkably, for Case I 
in Fig. BJc), a novel state occurs where at one sam- 
ple boundary spin-up and spin-down currents counter- 
propagate whereas at the other sample boundary spin- 
up and spin-down currents propagate in the same direc- 
tions. Using Landauer conductance formula for the four- 
terminal geometries we find that for Case I in Fig. [4^c) 
with a general N x , there is spin-polarized charge current 
(2N X + l/2)Ve 2 /h flowing from terminal 1 to terminal 
2 with a Hall voltage V/2, leading to a Hall conduc- 
tance a xy = (4JV X + l)e 2 /h. We note that this feature 
differs this state from the state of topological insulator 



since the Hall voltage in the latter state is 00, not V/2. 
For Case II in Fig. BJd) with a general N x , there is also 
a spin-polarized charge current (2N X + l)Ve 2 /h flowing 
from terminal 1 to terminal 2 but without a Hall voltage, 
leading to a resistance \/a xx = h/(2N x + l)e 2 . Inter- 
estingly, for Case I in Fig. SJc), there exist spin cur- 
rents (2N X — 3/2)Ve 2 /h flowing from terminal 1 and 
Ve 2 /h flowing from terminal 4, as well as a spin cur- 
rent (2N X — l/2)Ve 2 /h flowing to terminal 2, whereas for 
Case II in Fig.[3Jd), there exist a spin current 2N x Ve 2 /h 
flowing from terminal 2 to terminal 1, as well as a N x in- 
dependent spin current Ve 2 /h flowing from terminal 3 to 
terminal 4. We note that a wave-vector-dependent spin- 
filtering effect was also revealed recently by a calculation 
on the transport problem throu gh m agnetic barriers in 
graphene with Zeemann splitting [13J]. 

A natural question is that the SMF we considered here 
is ideal and changes direction abruptly, while in real con- 
ditions there exists a length scale I which is the distance 
covered by the SMF to change direction. Detailed cal- 
culation shows that if I is much less than the magnetic 
length Ib, our results will be independent of I. The / is 
estimated to be 1 ran, while for a typical magnetic field 
order of 10 T, Is ~ 10 nm, satisfying the criterion. 

In conclusion, graphene QHE under a modulated or- 
bital magnetic field has been investigated. The current- 
carrying edge states created by the modulated magnetic 
field give rise to a novel quantized Hall conductance, 
which can be checked by a standard four-terminal Hall 
measurement. By varying the UMF, the four-terminal 
graphene sample is expected to undergo a transition from 
a state with novel QHE to one without Hall effect. The 
effect of Zeeman spin splitting is also discussed and a 
novel state and its corresponding spin Hall currents are 
predicted. 
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